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O Abstract. We investigate the initial conditions problem for multifield inflation. In these 

-<-J scenarios the pre-inflationary dynamics can be chaotic, increasing the sensitivity of the onset 

C^ of inflation to the initial data even in the homogeneous limit. To analyze physically equivalent 

scenarios we compare initial conditions at fixed energy. This ensures that each trajectory is 

C\J counted once and only once, since the energy density decreases monotonically. We present 

K^ a full analysis of hybrid inflation that reveals a greater degree of long range order in the 

.^ set of "successful" initial conditions than was previously apparent. In addition, we explore 

CnI the effective smoothing scale for the fractal set of successful initial conditions induced by the 

'^ finite duration of the pre-inflationary phase. The role of the prior information used to specify 

■^ the initial data is discussed in terms of Bayesian sampling. 
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1 Introduction 

The standard hot big bang is synonymous with the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric. This model imposes maximally symmetric initial conditions on the metric 
and the mass-energy distribution, as specified on an arbitrary initial spatial hypersurface. 
These initial conditions acausally correlate spacelike-separated regions and require further 
fine-tuning for the Universe to be spatially flat at late times, leading to the well-known 
horizon and flatness problems. Famously, these problems are resolved by inflation [1-4], 
which grafts a phase of accelerated expansion onto the very early universe, setting the stage 
for the standard cosmology. During inflation the comoving Hubble volume contracts and the 
visible universe is driven toward the spatially flat FLRW universe. 

Given that inflation attempts to explain the otherwise ad hoc initial conditions of the 
standard hot big bang, a viable inflationary mechanism must itself be free of tunings. Tunings 
can appear as technically unnatural parameter values in the inflaton sector or the need for 
a special pre-inflationary field configuration: the latter question is the focus of this paper. 
Inflationary models with unnatural initial conditions are at best incomplete and, at worst, 
not viable as descriptions of the early universe. Moreover, the level of tuning required to 
ensure the onset of inflation can differ substantially between scenarios with largely degenerate 
observational predictions, providing a possible mechanism for discriminating between them. 

The initial conditions problem arises even in the purely homogeneous limit. For instance, 
chaotic inflation [5] begins for a large range of initial field values, but new inflation with 
the Coleman- Weinberg potential [3, 4] requires a special initial state. Inflationary models 
with multiple scalar degrees of freedom introduce a further level of complexity. With two 
or more fields the homogeneous dynamics are potentially chaotic, as first pointed out in 
Ref. [6] and also discussed by Refs [7, 8]. Chaos is synonymous with sensitive dependence 
on initial conditions, rendering multifield models qualitatively different from their single field 
counterparts. Multifield scenarios are widely studied and more natural in many settings. 
In particular, string theoretic inflationary scenarios often possess many scalar degrees of 
freedom. Further, even if a model has an effective single-field description once inflation is 
underway, the pre-inflationary phase may contain many interacting fields. Several analyses 
of the initial conditions problem for multifield inflation exist [7-14] and we return to this 
question here. 



Set Description 

X Initial conditions surface, however defined 

Z Set of initial conditions with zero velocity 

Ce Set of initial conditions with equal energy E 

Se Successfully inflating subset of Ce 

Te Non-inflating subset of Ce 

Be Boundary between Se and Te 

Table 1. Subsets of phase space: notation. 



We need to sample the "initial conditions space" X for these scenarios, determining the 
overall fraction that inflates and the topology of the inflationary region within this space. A 
homogeneous, spatially flat universe containing A*" scalar fields (^i with arbitrary interactions 
has 2N independent degrees of freedom since the scale factor can be eliminated by the 0-0 
Einstein equation. The solutions to the equations of motion — called trajectories or orbits 
— are non-intersecting curves that fill the 2A^-dimensional phase space. Yet, the initial field 
values and velocities are not independent or identically distributed (iid) random variables as 
different points in X are correlated by the solutions to the field-equations, i.e., many points 
belong to the same trajectory. (See Table 1 for a summary of our notation.) 

The phase space is foliated by surfaces of equal energy, Ce- The energy density p = E'^ 
of FLRW universes is monotonic, decreasing in a homogeneous universe as 
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where the Hubble parameter H cc E"^ and overdots denote derivatives with respect to coor- 
dinate time t. For a specific energy E, orbits intersect Ce once and only once, identifying 
each point on Ce with a unique solution to the equations of motion. To build a well-defined 
sample of trajectories we choose initial conditions from the constraint surface Ce- 

Many previous treatments of the multifield initial conditions problem [9-13] have been 
based on Z, the A^-dimensional subset of X on which all velocities vanish simultaneously. 
Although the set of infiationary trajectories intersecting Z is easier to sample than Ce, orbits 
for which all velocities vanish at the same instant are not generic, given the finite duration 
of the pre-infiationary era. Many orbits thus never intersect Z, while in principle others may 
intersect it multiple times, correlating apparently distinct points — issues that cannot arise 
when sampling from Ce- By contrast, Ref. [8] samples the full phase space X and also varies 
the parameters in the potential itself, effectively marginalizing over the energy scale E. The 
current paper is the first analysis of the initial conditions problem for multifield infiation that 
does not either (a) study a lower-dimensional surface in the initial conditions space (however 
defined) or (b) sample the entirety of X, disregarding the fact that different points belong to 
the same solution of the field equations.^ 

Any description of the primordial universe breaks down above some energy scale, and 
this scale defines the appropriate initial conditions hypersurface Ce for a well-specified 



^Tetradis [11] presents a single equal energy slice, although a degree of freedom was removed by requiring 
the velocities to be equal, similar to the projections we introduce for convenience in Figs 3-5. 



model. ^ This energy will be associated with some scale in the particle physics sector, such as 
the characteristic size of extra dimensions, the string scale, next-to-leading order corrections 
to Einstein gravity, or ultimately its breakdown at the Planck scale. Points on Ce are thus 
physically commensurate, whereas points in I span several orders of magnitude in energy. 
Qualitatively, we will also find that the set of successfully inflating points has a simpler struc- 
ture and more obvious long range order when chosen from Ce rather than Z, allowing us to 
better understand the underlying cosmological dynamics. 

Beyond the choice of initial conditions surface, we must also specify the prior probability 
distributions (in Bayesian terms) for the initial field values and velocities. If a; is a probability 
distribution that weights an initial condition xq according to how well its final state matches 
the observed universe, then the expected value of w over initial conditions xq G Ce is 

' ^(xo)P^(xo)d^xo^-^{u;(x«)} ,, , (1.2) 

Ce "- j=]^ ^ -^ ^0 ^-Pb 

where Pe is the prior probability distribution for initial conditions on the constraint surface 
Ce and the sum is evaluated at n points sampled from Pe- The prior in Eq. (1.2) acts 
as the probability density function for initial conditions on the space of FLRW universes. 
The form of Pe is only weakly constrained by fundamental considerations. The freedom to 
choose Pe is analogous to the measure problem in the multiverse [15-19], albeit restricted 
to the subspace of homogeneous FLRW universes. In many previous works the prior is often 
not directly discussed, and thus implicitly defined as a uniform distribution on the initial 
conditions. We consider several possible choices of prior (all of which are uninformative) 
and vary the energy of the surfaces Ce- We find that the choice of prior significantly alters 
the fraction of trajectories that lead to inflation, potentially distorting conclusions about the 
extent to which a given inflationary model requires fine-tuned initial conditions. 

In what follows we work with a widely studied two-field model: canonical hybrid or false- 
vacuum inflation [20-22]. We relate the initial conditions problem to that of determining the 
(fractal) topology and the geometry of the subset of points Se C Ce that successfully inflate, 
since this is independent of the choice of prior. Like Refs [8, 13] we see that Se has a fractal 
topology due to the presence of chaos in the underlying dynamical system, demonstrating 
that hybrid inflation has regions of phase space where orbits are highly sensitive to their 
initial conditions and confirming the results of Ref . [6] . Hybrid infiation is associated with a 
blue power spectrum^ at odds with recent astrophysical data [25-29]. However, our primary 
focus is not hybrid infiation itself, but developing tools that can be used to understand 
the initial conditions problems in generic models of multifield inflation. We use this model 
because (a) it is the prototypical multifield model with chaotic dynamics and a narrowly 
defined inflationary attractor; (b) we are primarily interested in the onset of inflation; and 
(c) to make contact with previous work. 

The pre-inflationary universe is dissipative, so the fractal structure must have a nontriv- 
ial scale dependence: there is necessarily a minimum scale below which two nearby trajecto- 
ries will remain correlated until they reach either the inflationary attractor or a minimum of 
the potential, smoothing Se below this scale. Conversely, while we assume classical homo- 
geneity, quantum fluctuations prevent the universe from being perfectly smooth. If Se has 



^If the potential has one or more local minima where Va > 0, choosing _E* < Va will necessarily exclude 
all trajectories which evolve toward these minima. 



^Although see Refs [23, 24]. 



structure on scales smaller than a typical fluctuation we cannot sensibly define the homoge- 
neous limit for this system. Consequently, we propose a sampling technique that identifies 
regions where Se has structure below this minimum scale. 

The paper is arranged as follows: in Section 2 we review hybrid inflation and discuss 
its dynamics. In Section 3 we describe our numerical methods, characterize the properties of 
the set of inflationary trajectories with different energies and priors Pe, and investigate the 
fractal dimension of Se- In Section 4 we discuss the implication of our results and identify 
future lines of enquiry. 

2 Inflationary dynamics 

For simplicity we consider two homogeneous scalar fields, ip and the inflaton <f), interacting 
through a potential V{ip,(j)) in a homogeneous FLRW universe. The equations of motion are 

■ dV ... dV 

(l) + 3H(j)+—- = and ^P + 3HiP + ^— = 0, (2.1) 

and the Hubble parameter H can be eliminated by the 0-0 Einstein equation 
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where Mpi is the Planck mass. Following Refs [6-13] we consider hybrid inflation [20-22] 
with the potential 
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with real parameters A, M, /i, and v. Inflation occurs in the "inflationary valley" with 
'(/' ~ and \(f)\ > (j)c, where (pc = v2ij'^/M is the critical point at which the effective mass 
of ip becomes complex and inflation comes to an end. The potential is symmetric under 
(f) — )■ —(j) and ip —7- —ip, with two equivalent valleys for (p > (pc and (f) < ~4'c and minima at 
{tp,4>} = {ibM, 0}. Orbits will either enter one of the false-vacuum inflationary valleys or 
evolve directly toward one of the true vacua. 

We set the amplitude Ag of the dimensionless power spectrum Vn to be roughly com- 
patible with the WMAP9 data [26, 27], which fixes the potential energy scale. This results 
in 

where ey = {M^j2){V^/Vf is the slow-roll parameter. Setting M = .03Mpi, // = 500 Mpi, 
and ly = .015 Mpi and assuming perturbations are generated when ip ^ and (p ^ (pc, we 
derive A ~ 6.8 x 10~^ -/Vfpi.^ Lastly, Ref. [30] determined that quantum fluctuations dominate 
the classical field evolution in the inflationary valley when 

A > Aq = 4^V3M|i^. (2.5) 

For our parameters A„ = 9.6 x 10~ Mpi, so the classical equations are self-consistent. 



The super-Planckian value of fi is an artifact of this definition of the potential; the actual mass term is 
ml = 2A*/ti^ ^ 10"^'' M|i, and safely sub-Planckian. 
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Figure 1. Distribution of successfully inflating initial conditions drawn from the zero- velocity slice 
Z. White areas have the highest number of successful points; darker regions have the fewest. This 
plot matches Fig. (1) of Ref. [8], with M = .03Mpi, /x = 636 Mpi, and v = .0173 Mpi. The figures 
are similar (verifying our codes and algorithms) but are not expected to be identical, due to different 
binning procedures. 



The inflationary valley is a small subset of the total phase space, which might suggest 
the model has a fine-tuning problem. References [9, 11, 12] considered sub-Planckian initial 
field values on the zero- velocity surface Z^ pessimistically concluding that — in the absence 
of effects that increase the friction experienced by the fields — only trajectories which start 
inside the inflationary valley yield 60 e-folds of inflation. By contrast, Ref. [10] was more op- 
timistic, showing that a supergravity-inspired hybrid inflation model has a significant number 
of "successful" points outside the inflationary valley. With more exhaustive sampling of Z^ 
subsequent studies by Clesse, Ringeval, and Rocher [8, 13] extended this optimistic conclusion 
to the potential (2.3). They showed that successful initial conditions are distributed in an 
intricate series of patches and fine lines outside the inflationary valley, with a fractal bound- 
ary separating inflating and non-inflating initial conditions. The distribution of successfully 
inflating initial conditions on Z^ for a specific scenario from Ref. [8], is reproduced in Fig. 1. 
The fine-tuning problem may also be less serious if the initial field values are assumed to be 
super-Planckian or if the interaction term dominates [13, 31]. 

References [8, 23] also present a Markov Chain Monte Carlo (MCMC) sampling of 
all possible parameter choices and sub-Planckian field configurations, including those with 



initial velocities. The conclusion was that 60 e-folds of inflation is generic for the potential 
in Eq. (2.3) and flne-tuned initial conditions in the inflationary valley are not required. 
Although we have argued that sampling from any two-dimensional subspace, such as Z^ is of 
limited benefit, sampling the whole four-dimensional space X may not be strictly necessary, 
even though an MCMC technique marginalizes the unknown initial energy E. We instead 
choose to explore how fine-tuned the initial conditions must be when sampling from constraint 
surfaces Ce that incorporate the energy constraint Eq. (1.1). 

3 Numerical results 

3.1 Method 

We numerically integrate Eqs (2.1)-(2.3) using a backward-difference formula implemented by 
the FcvODE package from the Sundials computing suite [32] . We sample initial conditions 
from constraint surfaces Ce with constant energy density 

P=^V'o + ^<Ao + nV'o,<Ao) = i^', (3.1) 

where 

^ = 10*Mpi for ie {-5, ...,0}. (3.2) 

The last 60 e-folds of infiation occur at {V','/'} ~ {0,(/)c} with E ~ 10"^Mpi. With E = 
lO'^Mpi we are at the limit of classical Einstein gravity; we only include this case to illustrate 
the underlying dynamical system. 

We stop integrating when either (a) the orbit achieves more than 60 e-folds during 
infiation or (b) p < b^ and the trajectory is trapped by the potential wells at {V')'/'} = 
{zbM, 0}. Initial conditions which lead to 60 e-folds of inflation are "successful" and define 
the subset Se-, while its complement — the "failed" points — comprise the subset Te^ The 
boundary between these sets, whose properties determine the extent to which they "mix," is 
denoted Be- 

We select points randomly on the constraint surface as follows. We first draw (^o and 
'(/'o from the uniform distribution over < {"^Oi^o} ^ •2Mpi, excluding any choices with 
y((/)o, V'o) > -E'^- The symmetry of the potential (2.3) allows the restriction to positive field 
values, whereas the upper bound is set to be consistent with Ref. [8]. Given these initial field 
values, the kinetic energy is typically dominant unless S ~ A ~ 10~^ M-p\. We apportion the 
remaining energy by drawing one of f i G {^O; V'o} from a uniform prior on the range 

- ^J2(E^ - Fo) <vx< ^2(^4 _ y^) (33) 

and giving the leftover energy to the other field velocity vi with the overall sign again chosen 
randomly. We maintain the symmetry between 0o and V'o by alternating the order in which 
these velocity terms are set.^ This procedure implicitly defines the initial priors Porig on the 
constraint surfaces Ce- 

To estimate the size of inhomogeneous fiuctuations at the initial energy E^ we note that 
b4> ~ H/2tt and H ~ E'^/Mpi for a massless field in de Sitter space. Similarly, the minimal 



^ Again, additional constraints can be added, e.g. data matching for rig, r, or other observables. 

^This prior generates the high tails in the velocity distributions seen in Fig. 6. If the first velocity chosen 
is «i, the second will be V2 = ±\/2{E'^ — Vo) — vf. Since vi is uniformly distributed, V2 is a quadratic 
distribution, favoring higher values. 
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Table 2. Total fraction of successfully inflating points sampled from priors Porig on the equal energy 
slices Ce — both excluding (/true) and including (/total) false positives from Se- Also shown are the 
number of successful points risucc, the number of false positives rifaiso, and the combined number of 
fail points and false positives nfaii. (f) The sampling procedure deviates from an "equal-area" sample 
as E ^ K. 

variation in velocities is expected to be of order H^ across a Hubble volume [30, 33, 34]. 
The fields ip and (j) ^'^^ iiot massless and the pre-inflationary universe is not de Sitter, but 
we can use this relationship to put an approximate lower bound on the homogeneity of the 
primordial universe. In regions of Se whose typical scale in any phase space dimension is 
less than 



A = {5iIjq,5(J)q,5i^jq,54>q] 
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the homogeneous approximation breaks down and further analysis is invalid or ambiguous. 
Physically, in these regions we cannot self-consistently assume that the primordial universe 
is homogeneous. 

We exclude these regions from Se by sampling Ce in clusters. We first choose points 
from Ce and integrate Eqs (2.1)-(2.3). At each point that successfully inflates we randomly 
draw 100 points within A of that point. ^ If any of these new points do not inflate, we conclude 
that the original point was a "false" (or perhaps ambiguous) positive. Although this approach 
incorporates points lying near (but not actually on) our designated equal energy surface, we 
do not weight our conclusions by these secondary points. 

This analysis does not address the inhomogeneous initial conditions problem; it sim- 
ply limits the extent to which the initial conditions can be self-consistently fine-tuned in a 
homogeneous universe, given that the chaotic dynamics of the potential may cause closely 
correlated trajectories to diverge exponentially. Fig. 2 shows three solutions of Eqs (2.1)- 
(2.3), at E = 10-^ Afpi with initial field values which differ by only 10~^Mpi. They eventually 
diverge, with each trajectory reaching a distinct end-state. If an inflation model has a frac- 
tal Se or Be that is distributed in a complex manner over Ce, then almost all successfully 
inflating initial conditions may be within A of an initial condition which does not inflate. 



3.2 Successful inflationary trajectories 

The fraction of successful points at any given energy E is summarized in Table 2, both 
including and excluding the "false positives." The highest probability for success is at higher 
energies. We should expect this since, given that the effective equation of state is not the 



^For the lowest values of E, < Sipo/jpn-iayc, 5(j!)o/0max, Sipo /ipmux, 50o/0max ^ ~ 10 ^'\ whicfi Is far below the 
resolution of our figures. We confirmed the accuracy of the Fcvode integrator in this domain by using an 
arbitrary precision integrator from Mathematica. 
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Figure 2. Parametric trajectory plots for three orbits at energy E = 10^^ Mpi initially separated 
by lO^^A/pi in the field values and 10~*Mp[ in the field velocities. These orbits are exponentially 
diverging, with each of the three trajectories branching at i sa 3600 A/pj^; one goes to each of the 
global minima a,t ip = ztM and one inflates. 



same on all trajectories, orbits accumulate on the narrow inflationary attractor over time, 
leaving a larger flux of orbits through the attractor at lower energies. In comparison to a 
sample drawn from Ce, an identical sample from a slice Ce' with E > E' will not place 
as much weight on trajectories inside the inflationary attractor and we expect to see fewer 
successfully inflating initial conditions on the lower energy surface. 

Figures 3 to 5 show two dimensional slices of Se at different values of E. We bin Se 
on a 1000 x 1000 grid: white regions are those with the highest number of successful points 
and the darkest regions have the fewest. We accumulated 2.5 million successful points on 
each slice and see minor stochastic variation in the number of points per bin. Figures 3, 4 
{E = 10~^Mpi), and b {E = 10~^ Mpi) show two-dimensional slices of Ce on which the 
initial velocities have equal magnitude \4)q\ = \iPq\- Looking at Figs 3 through 5 we can 
see areas where Se and J-e niix together, forming an intricate substructure similar to that 
seen in Refs [8, 13]. In Fig. 4 we also see contiguous regions and thick bands which reliably 
inflate and survive the subtraction of the false positives from Se- Qualitatively, Se exhibits 
considerable long range order when compared to Fig. 1. At higher energy, contiguous regions 
in Se occupy a larger portion of C^;, which can be seen clearly in Figs 3 and 4. The transition 
from Figs 3 to 5 shows how the geometry of Se changes with E. 

Figures 3 to 5 are projections of three dimensional regions and suppress information 
about the field velocities of the successful initial configurations. Intuitively, the points most 
likely to inflate (for given (J)q and V'o) would be those which had a large |i;^o| and small |'(/'o|- 
These points are essentially "launched" up the inflationary valley, while the slope of the 
potential focuses them toward smaller values of ^. To show this dependence on the initial 
values 00 and V'O) Fig. 6 shows histograms of these values sampled from the whole of Ce- The 
fraction of sampled points in Se as a function of initial velocity confirms our intuition: most 
successful points have larger \(J)q\ and smaller I'f/'ol- Points for which 0o ~ are particularly 
disfavored, further suggesting that the zero-velocity slice Z is unrepresentative of typical 
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Figure 3. Two dimensional slicings oiCE for E — 10 "^ Mp\, including the ambiguous, "false positive" 
points in Se- Parameters are A = 6.8 x 10"^ Mpi, M = .03Mpi, ^ = 500 Mpi, and v = .015 Mpj. The 
light and dark areas are regions that have a higher and lower density of points in Se, respectively. 
The results have been binned over a 1000 x 1000 grid. All velocities are of equal magnitude, however 
the left column has > 0, the right column is at </) < 0, the top row has i/' > 0, and the bottom row 
has ijj < 0. 



inflationary trajectories. We can also see the impact of the "false positives" in these plots: 
these are more frequent at high energies and in the limiting case E = Mpi all naively-inflating 
initial conditions are false positives, since A encompasses the whole of X in this limit. 

With sub-Planckian initial field values the kinetic energy dominates the potential energy 
for E ^ A. Thus, even if a trajectory starts inside the inflationary valley, its velocity is such 
that it is unlikely to remain there. For example, when E = 10~^ Mpi the slices of Se 
in Figs 3 and 4 show no special preference for points within the valley. In contrast, with 
E = W~^ Mpi ~ lOA the valley is clearly distinguishable, as shown in Fig. 5, but only when 
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Figure 4. Two dimensional slicings of Cg for E = lO^'^ A/pi, excluding any ambiguous or "false" pos- 
itives from the set of successfully inflating initial conditions, Se- All velocities are of equal magnitude; 
the left panel has -00 < and (j>o > 0; and the right panel has {tpo, 4>o} < 0. 
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Figure 5. Two dimensional slicings of Ce for E = 10~^ Mpi, where the checked background has 
l^(V'o, 4>o) > E'^ and has not been sampled. All velocities are of equal magnitude; the left panel has 
{■00, ^o} > and the right panel has {-00, 0o} < 0. 



the initial velocity of (j) is directed "uphill," i.e. with </>o > 0. Conversely, each slice contains 
many successful points that lie outside the inflationary valley. 

In Fig. 7 we project specific representative solutions of Eqs (2.1)-(2.3) onto the {V'l'A} 
plane for initial conditions with energies E = lO^^Mpi and E = 10~^Mpi. Trajectories 
which unambiguously inflate show little topological mixing and are all reflected off of the 
maximum of the potential T^nax = ^|i/)=0 toward the inflationary valley. For E > 10~^Mp 
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Figure 6. Histograms of (/)o (left column) and ■00 (right column). The rows have energy 10~^Mpi 
(top), 10^"^ Afpi (middle), and 10~^ Mpj (bottom). The gray background is the total sample from Ce 
and the blue foreground is the successful subset Se- M E = Mpi (not displayed), Se contains false 
positives only. 



most trajectories contain regions in which the field values are super-Planckian. We do not 
exclude these trajectories, but we could easily add this as a separate requirement for a viable 
inflationary scenario, in which case almost no successful inflationary trajectories exist at these 
energies. 

To quantify the sensitive dependence on initial conditions independently of our sampling 
procedure, we use the box-counting method to compute the fractal dimension of both Se and 
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Figure 7. Successfully inflating trajectories projected onto the {ip, cj)} plane, with initial conditions 
{^po,(j)o} marked in red. The left panel is at energy E — 10^^ A/pi and the right panel is at energy 
E ~ 10~* Mpi. The gray, checked region is where the magnitude of the field values exceeds Mp\. 
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Figure 8. The left panel shows the box counting of Be for E = 10^^ Mp\. The slope of the best-fit 
line d — 2.558 is the box counting dimension. The right panel is the box counting dimension d versus 
the energy E for sets Se and Be at the energies in Eq. (3.2). The red line with boxes indicates Se 
and the blue line with circles is Be- 

its boundary Be, including the "false" positives [35]. We first cover Ce with progressively 
smaller four- dimensional boxes of size 6, then count the number N{6) of (5-sized boxes in each 
covering. The box-counting dimension 



d 



5^0 log(l/(5) 



(3.5) 



is estimated by the slope of the line fitted to the linear portion of the curve log(A^) as a 
function of log{l/5). To compute the dimension of Be we count boxes that contain elements 
of both Se and J-e- Figure 8 shows both a typical fit (for Be with E = 10^^ Mpi) and the 
computed values of d for Se and Be- The result is sensitive to the detailed fitting procedure, 
which we trained by testing the algorithm on sets with known dimension, such as Cantor 
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-forig 


-1 square 


P, 


P^ 


E [Mpi] 


/true 


/total 


/true 


/total 


/true 


/total 


/true 


/total 


10° 


0.000 


0.498 


0.000 


0.497 


0.000 


0.704 


0.000 


0.294 


10-1 


0.447 


0.479 


0.426 


0.471 


0.651 


0.681 


0.243 


0.276 


10-2 


0.389 


0.434 


0.347 


0.408 


0.582 


0.624 


0.196 


0.242 


10-3 


0.300 


0.331 


0.245 


0.285 


0.438 


0.469 


0.162 


0.194 


10"^ 


0.205 


0.211 


0.188 


0.195 


0.247 


0.253 


0.163 


0.169 


10-5 


0.225 


0.225 


0.215 


0.215 


0.265 


0.265 


0.183 


0.183 



Table 3. Fraction / of sampled points from Ce that inflate — both excluding (true) and including 
(total) false positives. The sampling techniques (/'orig, /^square, Px, and P; ) are explained in the text. 



dust. Moreover, the regions considered here are multifractal, in that the dimension of 5^; will 
be a function of both position in X and the overall scale. The first is easy to see: each surface 
Ce contains regions in which essentially all points inflate (in these regions d ~ 3) and regions 
that are approximately isolated points (with d < 3). Consequently, the computed value of 
d is effectively a weighted average of at least two different regions, which explains why the 
dimension of Se is close to 3 but still measurably non- integer. Secondly, at very small scales 
Se must consist of smooth contiguous regions and on these scales we expect d — )• 3. These 
regions exist in spite of the chaotic dynamics due to the dissipative terms in Eqs (2.1)-(2.3) 
and put a lower limit on the mixing scale. 

3.3 The role of the prior 

It is well-known [15-19, 31, 36] that probability measures on different hypersurfaces result in 
different conclusions regarding the likelihood of inflation; we explore here how this relates to 
the choice of sampling prior. Although surfaces with different energies (as well as different 
initial conditions surfaces, such as a slice of constant comoving time) are homeomorphic to 
Ce and, by definition, have the same topology, the prior on Ce is not a topological property 
and is not preserved under either homeomorphism or a change of variables.^ Each initial 
condition surface then has a different prior and different likelihood for inflation, even given 
the same sampling technique. 

In Table 3 we compare different uninformative priors, defined implicitly through four 
sampling algorithms, on surfaces Ce at the energies in Eq. (3.2). Since the kinetic en- 
ergy is initially dominant for the energies and ranges we consider, we leave the selection 
method for {ipQ^cpo} the same as in Section 3.1, but vary the way we set the velocities 
V S {V'Oj'/'o}- The original prior Porig draws one velocity vi from a uniform distribution, 
bounded by ibi/2(£'^ — Vq), and then sets V2 by the energy constraint (3.1). All signs are 
chosen randomly and this procedure is alternated on subsequent choices of points to obtain a 
symmetric distribution in the velocities. The second prior Pgquare is similar, except we draw 
the square of the velocity vf from a uniform distribution bounded below by zero and above 
by 2(£'^ — Vo), with the sign of vi chosen randomly. Again, we alternate this to obtain a 
symmetric distribution. With this modest change in prior, the calculated fraction /true of Ce 



A homeomorphism is provided by time-translation along the integral curves of the equations of motion. 
Note that Z, being of a lower dimension, is not homeomorphic to Ce- 
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that inflates (excluding false positives) differs by only a few percent, with Psquare giving a 

slightly lower fraction at each energy. The fraction /true again decreases with decreasing E. 

We compare this to two priors P; and P: that are asymmetric in the velocities. For 

P; we always draw tpo from a uniform distribution bounded by ib-^2(£^^ — Vq) and always 

set 00 by the energy constraint. For P; we do the opposite: draw (j)o and set ipQ. This 
gives a uniform distribution in the sampled velocity vi , but a high-tail distribution similar to 
Fig. 6 in the velocity V2- The prior P; focuses more of the sample around ipQ ~ 0, the area 

identified as being most-likely to inflate, and P; gives more points around </>o ~ 0, the area 
least-likely to inflate. At -E = Mpi the difference in /total (including false positives) between 
the asymmetric priors is as much as 41.0 percentage points. The differences decrease with 
decreasing E, indicating that later-time hypersurfaces become progressively independent of 
the prior. However, Table 3 demonstrates how any measure of / is prior-dependent, especially 
with respect to the implicit dependence of the prior on the initial energy. 



4 Conclusion 

We have considered the initial conditions problem for multifield inflation, quantifying the 
likelihood of inflation by sampling an initial conditions surface, evolving the points numeri- 
cally, and dividing them into successfully and unsuccessfully inflating sets. We draw initial 
conditions from an equal energy slice of phase space, denoted Ce, the maximum energy at 
which the underlying theory is assumed to be an accurate description of the primordial uni- 
verse. Since FLRW universes have a monotonic energy density, sampling initial conditions 
from Ce ensures that we count only unique solutions to the equations of motion. A sample of 
points from Ce is thus a well-deflned sample of homogeneous universes. Typically, we cannot 
predict the flux of orbits through Ce and must choose a prior, accordingly. We considered 
four different uninformative priors on Ce and showed that the likelihood of inflation varies 
by as much as a factor of roughly two between candidate priors. However, one can imagine 
scenarios where the prior dependence was much more dramatic. 

After specializing to hybrid inflation we examined the topology of the set of successful 
points Se, which is independent of continuous deformations to the prior. We confirm that 
both Se and the boundary between the successful and unsuccessful points is fractal for all 
sampled energies. The structure of Se, as seen in Figs 3 to 5, is qualitatively smoother 
than when initial conditions are chosen from the zero- velocity slice shown in Fig. 1. Further, 
since the equations of motion (2.1)-(2.3) are dissipative, there must be a small-scale cutoff 
to any fractal structure. However, quantum fiuctuations put a fundamental lower limit on 
the homogeneity of the early universe: if Se has structure below this scale, the assumption 
of homogeneity is not self-consistent. Fluctuations are larger at higher energies and above 
some critical energy E the number of viable, homogeneous scenarios is vanishingly small, even 
though the naive counting statistic suggests that a nontrivial fraction of the initial conditions 
space is infiationary. 

Our specific calculations are performed for the hybrid potential (2.3), but our underlying 
goal is to develop tools that can be applied to the initial conditions problem associated 
with generic multifield scenarios. Recent progress has been made by studying both random 
multifield models [37-42] and inflection point models [36, 42-46]. These approaches yield 
contrasting conclusions regarding the distribution of inflationary trajectories; applying the 
methods developed here to these models will be an interesting extension of this work. 
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This analysis assumes that the universe is initiahy spatiahy flat and homogeneous, but 
even if inflation begins without tuning in the homogeneous hmit there is no guarantee that 
this result will survive the addition of pre-inflationary inhomogeneities. Inhomogeneous pre- 
inflationary conflgurations were examined by Goldwirth and Piran [47-49], who showed that 
single-fleld chaotic inflation and new inflation [3, 4] remain robust in the presence of nontrivial 
inhomogeneity, provided that the initial fleld value is approximately correlated over several 
Hubble radii. We plan to examine this question for multifleld inflation in future work. 
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